Abrupt Changes in the Dynamics of Quantum Disentanglement 
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Entanglement evolution in high dimensional bipartite systems under dissipation is studied. Dis- 
continuities for the time derivative of the lower bound of entanglement of formation is found depend- 
ing on the initial conditions for entangled states. This abrupt changes along the evolution appears 
as precursors of entanglement sudden death. 
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Entanglement is a cornerstone of modern quantum 
physics [J, The evolution of entanglement in open 
quantum systems is a matter of increasing interest and 
new phenomena have been predicted [1, US B S B H 
[Tol | . One of the most outstanding effects arises when en- 
tanglement vanishes long before coherence is lost. It has 
been pointed out that systems composed of two qubits 
in a noisy environment can lose its entanglement in fi- 
nite time, a phenomena named as Entanglement Sud- 
den Death (ESD), even though full decoherence happens 
asymptotically. This feature appears for certain classes 
of states of two qubits under the action of independent 
reservoirs. Examples of these classes are the so-called 
"X" mixed states as well as some particular types of non- 
maximally entangled pure states [3]. 

The purpose of this work is to explore dynamical be- 
havior of entangled states in larger bipartite systems un- 
der the action of independent reservoirs. We show that 
unlike the case of two qubits, 3 <g> 3 systems may present 
not only ESD, but also intermediate abrupt changes in 
the disentanglement dynamics, i.e. the rate in which a 
given state loses its entanglement may change through- 



out the dissipative process even though coherence is lost 
in a constant rate. We show that these rate changes are 
associated with sudden changes in the rank of the par- 
tially transposed density matrix, which also provides an 
explanation for the sudden death of entanglement. 

We analyze the disentanglement dynamics for differ- 
ent initial states and show that abrupt changes may be 
present or not depending on the variation of a small num- 
ber of parameters. We also recover the result for two 
qubits when preparing the initial state in a 2cg> 2 subspace 
of the whole system. Finally, we interpret these results in 
terms of changes in the set of entanglement witnesses ap- 
propriated for the characterization of the entangled state 
in each part of the dynamics. 

In this work we use a general measurement for the 
lower bound of Entanglement of Formation (EOF) for a 
mixed state in m®n dimensions, which has been recently 
proposed [Til ] . This proposal is based on the comparison 
between two major criteria: (i) thepositivity under par- 
tial transposition (PPT criterion) [H,[l3| and (ii) the re- 
alignment criterion [3, EH ■ EOF for m g> n-dimensional 
systems (m < n) is defined as [TT|] 
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where m is the dimension of the first subsystem and 7 is 
given by 
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with A = max(||/9 T - 4 || , ||i?(p)||) and ^(x) = — x\og(x) — 
(1 — x) log(l — x). The trace norm II • II is defined by \\G\\ = 
tr(GG')i . The matrix p TA is the partial transpose with 
respect to the subsystem A, that is, pj^ji = Pjk,ti, and 
the matrix R{p) is defined as R{p)ij,kl — Pik,jl- 

The PPT criterion says that p TA > for a separable 
state [Hj]. On the other hand, the realignment criterion 



says that a realigned version of p, for a separable state 
must satisfy the condition: ||-R(/o)|| < 1. These conditions 
state that entanglement exists for A > 1. The maximum 
values that A(t) can assume depend on the dimensions of 
the bipartite systems. For example, for a maximal two 
qutrits entangled state A = 3, for two qubits A = 2. The 
minimum value for a separable state is always A = 1. In 
this work we use this quantity to study the time evolution 
of entanglement in the presence of dissipation. 

Let us consider entangled quantum states of two 
qutrits, with at most two excitations, in the presence 
of dissipation at zero temperature. Such situation can 
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FIG. 1: Evolution of A(t) for the initial mixed state of Eq.( 
[4]) with A = 0.1 as a function of the adimensional time Ft. 



be conveniently described by the evolution equation: 



P = T 2cipc * ~ c * 
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where Cj , c] describes annihilation and creation operators 
for bosonic modes and p is a 3 (8) 3 density matrix in the 
computational basis {| 0), | 1), | 2)} ® {| 0), | 1), | 2)} of 
both qutrits. 

Let us consider at first glance a class of initially mixed 
states of two qutrits which corresponds to a modification 
of a maximally entangled state given as follows: 
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where A is a real parameter ranging from < A < 1. 
In the extreme cases, A = 0, we have a separable state 
whereas for A = 1 we have a maximally entangled state. 
The Eq.Q can be solved for arbitrary decay constants, 
but for simplicity we reduce the problem to the simplest 
case Ti = T2 = T. By a numerical calculation we realize 
that ||p Ta || > f° r an times, so that, we need to 

concentrate only in A(i) = \\p TA \ \- 

Fig. [1] shows the evolution of A(t) for A = 0.1. As we 
observe, A(i) undergoes sudden changes along its evo- 
lution exhibiting discontinuous derivatives, finally evolv- 
ing to a situation where entanglement abruptly dies. As 
compared with the case of two qubits a reacher dynamical 
behavior of entanglement appears. From the definition of 
A(t), this feature must be closely connected with the tem- 



poral dependence of the eigenvalues of M = p 
In our case both analytical and numerical calculations 



of the eigenvalues of the matrix M can be carried out. 
From numerical calculations we realize that the abrupt 
changes of entanglement evolution are dominated by the 
behavior of a restricted number of eigenvalues given by: 
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are the density matrix elements pijjd- These eigenval- 
ues ([5]) are plotted in Fig. [2j where we observe that the 
times where they vanish are in exact agreement with the 
times where abrupt changes in the entanglement evolu- 
tion appear. From Eqs. (O, these times can be analyti- 
cally calculated in terms of the parameter A: 

t x = ln(2/(2-A)),t 2 = In(l/(1-A)),t 3 = ln(l/(l- VA)). 

(7) 

Fig. [3] shows the smooth behaviors of these times as 
a function of the parameter A defining particular two 
qutrits mixed states. From this picture we realize that 
the abrupt changes in the dynamics of entanglement will 
appear for all values of A in the interval [0,1]. In par- 
ticular for the maximally entangled state with A = 1, 
there is sudden change for t\ — In 2, and the time of the 
second and third sudden change, which is the ESD, goes 
to infinite, showing that the entanglement decays asymp- 
totically. Note that this result differs substantially from 
its two-qubit counterpart where the corresponding max- 
imally entangled states disentangle smoothly [7}. Also 
note that these abrupt changes can be mathematically 
interpreted as discontinuities of the derivative for the ex- 
pression A(t) = Y^i=i V Ei (*) > an d a sudden change in 
the evolution of A occurs whenever one of the nine eigen- 
values Ei becomes zero, as observed in Fig. [5] 

The analysis to explain these particular sudden 
changes in the evolution of entanglement has been done 
in terms of the eigenvalues of the matrix M. However, 
we can also understand these abrupt changes in A(t) by 
observing the behavior of the eigenvalues of the partial 
transpose matrix p TA . In our case only three eigenval- 
ues give us information about these sudden changes and 
are plotted in Fig. 2] We notice that these eigenvalues 
change from negative to positive values for specific times 
which are in agreement with the sudden changes in the 



0.0 0.1 0.2 0.3 0.4 0.5 

n 

FIG. 2: Relevant eigenvalues of matrix M for the case shown 
in Fig. Q] 
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FIG. 3: Times for sudden changes in the dynamics of entan- 
glement as a function of the parameter A. Solid line corre- 
sponds to the first sudden change, dashed line corresponds to 
the second sudden change, and dot-dashed line corresponds 
to the sudden death. 

entanglement evolution. In other words, the disentan- 
glement rate changes whenever the rank of the partially 
transposed matrix changes abruptly. We can also asso- 
ciate to each eigenvalue of p TA a corresponding entan- 
glement witness operator such that a, (t) = Tr (Wtp (t)) 
with i = 1,2,3 and each Wi is given by 

Wi = ± [|21)(21| - |11)(22| - |22)(11| + |12)(12|] , 
W 2 = | [|10)(10| - |00)(11| - |11)(00| + |01)(01|] , (8) 
W 3 = \ [|02)(02| - |00)(22| - |22)(00| + |20)(20|] . 

At t = 0, all three operators can be used to identify 
entanglement in p. As time goes by, they consecutively 
lose this capacity until there is no entanglement left. This 
suggests a geometrical interpretation to the phenomena 
here described which will be explored in further publica- 
tions. 

It is interesting to compare the case analyzed previ- 
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FIG. 5: Entanglement evolution for the case shown in Fig 
JT]) (dashed line) for A = 0.15 and for a initial state in the 
subspace {| 11), | 22)} (solid line) for x = 0-2 

ously with that of an initial state restricted to a two- 
dimensional subspace p = (1/2) (|11> (11 1 + |22)(22| + 
X|H> (22| +x|22)(ll|). Fig.dlshows the evolution of the 
entanglement for \ ~ 0.2 as compared with the state 
in Eq. 2] for A = 0.15. We observe that in the case of 
the initial condition restricted to two dimensional sub- 
space we have only one abrupt change in the evolution 
corresponding to the ESD which resembles the behavior 
observed for two qubits 0] . If we look at the eigenvalues 
of M we see that at the time when A goes to zero there 
is also an eigenvalue that goes to zero, indicating that an 
abrupt change occurs. A similar conclusion could actu- 
ally be obtained when looking at the negativity instead 
of Concurrence for the case of two qubits. 

In addition we could explore the entanglement evolu- 
tion for initial non maximally pure entangled states, for 
example, | $) = a | 00) + /3 | 11) + 7 | 22). For the 
sake of simplicity we consider a, (3 and 7 real positive 
numbers. In this case, a richer dynamics for the entan- 
glement can be observed. Depending on the choice of the 
amplitudes we can have asymptotic decay, sudden death, 
sudden changes or a combination of them. Times corre- 
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From these expressions we realize that the entanglement 
dynamics exhibit: (a) asymptotic decay for (a > /3 > 7), 
(6) one sudden change and asymptotic decay for (/3 > 
a > 7, or a > 7 > /3), (c) two sudden changes and 
asymptotic decay for {(3 > 7 > a), and (d) two sudden 
changes and ESD for (7 > (5 > a). Fig. [6] shows two 
particular dynamics evolution for the cases (b) and (d). 
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FIG. 6: Entanglement evolution for the non maximal pure 
state. Dashed line corresponds to a = 0.2386, P = 0.9545, 
7 = 0.1790 and solid line corresponds to a = 0.1790, /3 = 
0.2386, 7 = 0.9545. 



sponding to the sudden changes and the ESD time are 
given by: 



ti = -£ln(l-£),t3 = -£ln(l-f), 
t 2 = -fin 
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where 
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In summary we have studied the evolution of entangle- 
ment for high dimensional dissipative quantum systems. 
By evaluating the entanglement contained in the system 
using the Chen, Albeverio and Fei measure we have ob- 
served outstanding new effects. Quantum correlations 
undergo abrupt changes as precursors of ESD. These can 
be characterized by observing the eigenvalues of the M 
matrix which defines the amount of entanglement for the 
quantum system. The dynamical changes are related to 
sudden changes in the rank of the Matrix M. Similar 
behavior can be found for both initially mixed or pure 
states and the ESD is recovered as a particular case of 
these sudden dynamical changes. 
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